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ABSTRACT 

An E 7 jet oc E' p ' relationship with a small scatter for current 7-ray burst (GRB) data was recently reported, 
where i? 7 ,jet is the beaming-corrected 7-ray energy and E' p is the vF u peak energy in the local observer frame. 
By considering this relationship for a sample of 12 GRBs with known redshift, peak energy, and break time of 
afterglow light curves, we constrain the mass density of the universe and the nature of dark energy. We find that 
the mass density £Im = 0.35±q [5 (at the lcr confident level) for a flat universe with a cosmological constant, and 
the w parameter of an assumed static dark-energy equation of state w = — 0.84±Qgj (la). Our results are consistent 
with those from type la supernovae. A larger sample established by the upcoming Swift satellite is expected to 
provide further constraints. 

Subject headings: cosmological parameters — cosmology: observations — gamma-rays: bursts 



1. INTRODUCTION 

Type la supernovae (SNe la) have been playing an important 
role in modern cosmology. Early observations of SNe la at red- 
shift z < 1 strongly suggest that the expansion of the universe at 
the present time is accelerating (Riess et al. 1998; Perlmutter et 
al. 1999). Since then, the nature of dark energy (with negative 
pressure) that drives cosmic acceleration has been one of the 
greatest mysteries in modern cosmology (for reviews see Pee- 
bles & Ratra 2003; Padmanabhan 2003). Recent observations 
of 16 higher-redshift (up to z — 1.7) SNe la present conclusive 
evidence that the universe had once been decelerating (Riess et 
al. 2004). These newly-discovered objects, together with pre- 
vious reported SNe la, have been used to provide further con- 
straints on both the expansion history of the universe and the 
equation of state (EOS) of a dark energy component (Riess et 
al. 2004; Wang & Tegmark 2004; Daly & Djorgovski 2004; 
Feng, Wang & Zhang 2004). 

7-ray bursts (GRBs) are the brightest electromagnetic explo- 
sions in the universe. It has been widely believed that they 
should be detectable out to very high redshifts (Lamb & Re- 
ichart 2000; Ciardi & Loeb 2000; Bromm & Loeb 2002). 7- 
ray photons with energy from tens of keV to MeV, if produced 
at high redshifts, suffer from no extinction before they are de- 
tected. These advantages would make GRBs an attractive probe 
of the universe. Schaefer (2003) derived the luminosity dis- 
tances of 9 GRBs with known redshifts by using two luminosity 
indicators (the spectral lag and the variability). He obtained the 
first GRB Hubble diagram with the la constraint on the mass 
density SIm < 0.35. 

A correlation between the isotropic-equivalent 7-ray energy 
(Eyjsa) and the vF v peak energy (E' p ) in the local observer 

1 /2 

frame, E p oc E^ iso , was discovered by BeppoSAX observations 
(Amati et al. 2002; Yonetoku et al. 2004), and confirmed 
by HETE-2 observations (Sakamoto et al. 2004; Lamb et al. 
2004). It not only holds among BATSE GRBs (Lloyd-Ronning 
& Ramirez-Ruiz 2002) but also within one GRB (Liang, Dai 
& Wu 2004). Its possible explanations include the synchrotron 
mechanism in relativistic shocks (Zhang & Meszaros 2002; Dai 



& Lu 2002) and the emission from off-axis relativistic jets (Ya- 
mazaki, Ioka & Nakamura 2004; Eichler & Levinson 2004). 
However, the dispersion around this correlation is too large to 
obtain useful information on the universe from the current GRB 
sample. 

Ghirlanda, Ghisellini & Lazzati (2004) recently found a 
new relationship between the beaming-corrected 7-ray energy 
(E-y jet) and the local-observer peak energy, £ 7 j e t oc E' 1 ' 5 , with 
a small scatter for current GRB data, suggesting that GRBs are 

a promising probe of the universe. In principle, this relation- 
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ship can be derived from the E' p oc E^ is0 correlation combined 
with the afterglow jet model. In this Letter, we constrain the 
mass density of the universe and the nature of dark energy by 
considering this relationship with a sample of 12 GRBs with 
known redshift, peak energy, and break time of afterglow light 
curves. We show that GRBs appear to provide an independent 
and interesting probe of fundamental quantities of the universe. 

2. SAMPLE SELECTION AND STANDARD CANDLES 

By searching for GRBs in the literature, we have found to- 
tal 14 bursts of which redshift z, observed peak energy E p , and 
break time fj of afterglow light curves are available. Table 1 
lists a sample of 12 GRBs, but the other two events, GRBs 
990510 and 030226, are not included. The reason is as fol- 
lows: the analysis of this Letter and Ghirlanda et al. (2004) is 
based on the afterglow jet model (Rhoads 1999; Sari, Piran & 
Halpern 1999). In this model, a relativistic jet, after emitting 
a fraction r] 7 of its kinetic energy into prompt 7 rays, expands 
in a homogeneous medium with number density of n. As the 
jet sweeps up more and more medium matter, its Lorentz fac- 
tor declines. When the Lorentz factor equals the inverse of the 
jet's half opening angle 9, the afterglow light curve presents a 
break. However, this model cannot well fit the afterglow data of 
these two bursts, because the predicted break spans about two 
orders of magnitude in time when light travel time effects are 
taken into account, and thus the theoretical light curve is too 
smooth to be consistent with the observed sharpness (Rhoads 
& Fruchter 2001; Wei & Lu 2002). For GRB 990510, the re- 
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quired spectral index of the electrons is less than 2, being incon- 
sistent with the shock acceleration theory (Wei & Lu 2002). In 
addition, the afterglow data of GRB 030226 suggest that its en- 
vironment might be a low-density wind rather than a constant- 
density medium, also conflicting with the model (Dai & Wu 
2003). 

According to the afterglow jet model (Sari et al. 1999), 
the jet's half opening angle is given by 9 = 0.161(1 + 

zT 3 ^l^E~yi 52 n /S r]\ /S , where E lMo . 52 = £ 7 ,i S o/10 52 ergs, 
tjj = tj/l day, « = n/\ cm" 3 , and ^ = 0.2 (Frail et al. 2001). 
Only for few bursts in Table 1 the medium density was ob- 
tained from broadband modelling of the afterglow emission 
(e.g., Panaitescu & Kumar 2002). For those bursts with un- 
known n, we assume the median density n ~ 3 cm" 3 as in 
Ghirlanda et al. (2004). The isotropic-equivalent 7-ray energy 
of a GRB is calculated by 
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where 5 7 is the fluence (in units of erg cm" 2 ) received in some 
observed bandpass and k is the factor that corrects the ob- 
served fluence to the standard rest-frame bandpass (1-10 4 keV) 
(Bloom, Frail & Sari 2001). For a Friedmann-Robertson- 
Walker (FRW) cosmology with mass density Q M and vacuum 
energy density Q\, the luminosity distance in equation (1) is 

d L = c(l+z)Ho 1 \n k \-V 2 sirm{\n k \V 2 

x / ^[(l+z) 2 (l + MZ )-z(2+z)0 A ]- 1 / 2 }, (2) 
Jo 

where c is the speed of light and Hq = 100/zkms" 1 Mpc" 1 is 
the present Hubble constant (Carroll, Press & Turner 1992). In 
equation (2), Q, k = 1 - Qm~ Qa, and "sinn" is sinh for > 
and sin for < 0. For = 0, equation (2) turns out to be 
c(l +z)Hq 1 times the integral. In this section, we assume a flat 
universe (i.e., = 0) because of both an expected consequence 
of inflation and the observed characteristic angular size scale of 
the cosmic microwave background fluctuations (Spergel et al. 
2003 and references therein). 

From equations (1) and (2), we obtain the beaming-corrected 
7-ray energy E 7 j et = (l-cos6)E Jt i SO , that is, 



-7 jet 



1.30xl0 50 ( l + z)- 3 /^XSo,5 2 ^ /4 < /4 ergs. 



(3) 



Figure 1 plots Ejj et versus E' p for the GRB sample listed in 
Table 1, with Q M = 0.27, Q A = 0.73 and h = 0.71. We find 
that 2? 7 j e t and E' are strongly correlated with a correlation co- 
efficient r s = 0.99 ±0.08 (with a probability of < 10" 4 ). The 
best fit is (£7 jet/ 1 50 ergs) = (1 . 12 ± 0. l2)(E' p / 100 ke V) 1 ■ 50±0 08 
with a reduced Xdof = 0.53. We note this power to be insensi- 
tive to Qm- In addition, although the peak energy E' p and the 
low-energy spectral index a in Table 1 appear to evolve with 
redshift (Amati et al. 2002), this evolution doesn't affect the 
above relation as shown in Figure 1 . These results imply that 
GRBs are standard candles. 

3. HUBBLE DIAGRAM AND COSMOLOGICAL CONSTRAINTS 

We first derive the observed luminosity distance from the 
GRB sample. Considering a relationship (E~, j et /10 50 ergs) = 
C(£'/100keV) L5 (where C is a dimensionless parameter), we 



d L = 2.37 x 10 : 



(l+z) 2 C 2 / 3 E ; , 
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cm, 



(4) 



where E p = E' p /{\ +z) is in units of keV. Thus, the observed dis- 
tance modulus of a GRB is ^ = 51og(<f L /10pc) with an error 
of 



= 2.17 



^A 2 + (^2.) 2 + (^) 2 + (^ 

E p J \2S 7 J \2tjJ \6n 
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where ue p , crs, t , a, , and a n are the errors in the peak energy, flu- 
ence, break time and medium density of the GRB, respectively. 

We plot a Hubble diagram of our GRB sample in Figure 2 
based on equations (4) and (5). This figure also presents a Hub- 
ble diagram of the current SNe la sample. Both Hubble dia- 
grams are consistent with each other. However, GRBs and SNe 
la have mean uncertainties of 0.09 and 0.05 in the log of the 
distance, respectively, and thus GRBs are about twice worse in 
accuracy than SNe. 

For an FRW cosmology with VIm and Qa, equation (2) gives 
the theoretical distance modulus ^th = 5 log^/lOpc). The like- 
lihood for these cosmological parameters can be determined 
from a x 2 statistic, where 



[Mthfc;/z,^M,^A)-Mob,i(Q] 2 



(6) 



We consider all possible values of the parameters h and C to be 
h e (0.68,0.75) (Bennett et al. 2003) and C £ (1.00, 1.24) (see 
§2). The confidence regions in the Qm-^a plane can be found 
through marginalizing the likelihood functions over h and C 
(i.e., integrating the probability density P oc e~ x I 1 for all val- 
ues of h and C). We plot contours of likelihood (from 1 a to 3cr) 
for unknown curvature fl k in Figure 3. As shown for a flat uni- 
verse, with the current sample, CIm < 0.62 (at the 2a confidence 
level), and the la contour contains the (Qm, ^a) = (0.27,0.73) 
point corresponding to the "concordance" model. We measure 
Q M = 035±° \ 5 5 (la). 

There are several alternative approaches to calculate the lu- 
minosity distance (also see Riess et al. 2004). We here consider 
a flat universe and a constant EOS, w = Pde / Pdec 2 , of a dark en- 
ergy component (Garnavich et al. 1998; Perlmutter et al. 1999). 
In this case we have 
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d L =d\+z)iv / dz[n M (i+zf 
+(i-n M )(i+z) MUw, r 1/2 - 



(7) 



Figure 4 presents contours of likelihood in the Qm~w plane (af- 
ter marginalizing over h and C). The solid contours consider 
a prior of VIm = 0.27 ± 0.04 by assuming its Gaussian distribu- 
tion, similar to Riess et al. (2004). We see w = -0.84±g;g (la), 
which is consistent with the value of w expected for a cosmo- 
logical constant (i.e., w = -1). 



4. CONCLUSIONS 



The^jetoc^ 
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relationship with a small dispersion was re- 
ported by Ghirlanda et al. (2004) and confirmed in this Letter. 
The advantages of considering this relationship as a probe of the 
universe are (1) that GRBs have been detected at redshifts up to 
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z ~ 4.5, (2) that 7 rays suffer from no extinction, and (3) that 
we don't think of luminosity evolution. These advantages led 
us to constrain the mass density of the universe and the nature 
of dark energy. We found that the mass density M = 0.35±[J [5 
(la) for a flat universe with a cosmological constant, and the w 
parameter of the dark-energy EOS w = -0.84±[] 33 (la). Riess 
et al. (2004) measured Vt M = 0.29±g;g| and w = -1 .02±°{| (la) 
for the current SNe la sample. Therefore, our results are con- 
sistent with those from SNe la. 

The upcoming Swift satellite with an energy range of 0.2- 
150 keV will be scheduled for launch in 2004 September 
(Gehrels et al. 2004). Swift is expected (1) to detect more than 
100 bursts per year, (2) to observe X-ray and UV/optical after- 
glows at times of 1 minute to several days after the burst, and 



(3) to detect very-high-redshift GRBs. Thus, it is expected that 
a larger sample of GRBs established by Swift provides a further 
probe of the universe. Such a probe opens up a new window on 
the cosmic distance scale far beyond the reach of SNe la. We 
call this research field GRB cosmology, corresponding to the 
well-known supernova cosmology. 
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Natural Science Foundation of Yunnan (2001 A0025Q), and the 
Research Foundation of Guangxi University. 
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Table 1 
Sample of 7-ray bursts 



GRB 


redshift 


E p (o Ep f 


a a 


/3 a 




range b 


fj(^) c 


n(o„Y 


ref 2 






(keV) 






(erg cm" 2 ) 


(keV) 


(days) 


(cm" 3 ) 




970828 


0.957 


297.9(59.3) 


-0.7 


-2.07 


9.6E-5(0.9) 


20-2000 


2.2(0.4) 


3(0.33) 


1,13,19 


980703 


0.966 


255.3(50.9) 


-1.31 


-2.39 


2.3E-5(0.2) 


20-2000 


3.4(0.5) 


28(10) 


2,13,20,20 


990123 


1.6 


780.8(61.9) 


-0.89 


-2.45 


3.0E-4(0.4) 


40-700 


2.04(0.46) 


3(0.33) 


3,14,21 


990705 


0.843 


188.8(15.2) 


-1.05 


-2.2 


7.5E-5(0.8) 


40-700 


1.0(0.2) 


3(0.33) 


4,14,22 


990712 


0.43 


65.0(10.5) 


-1.88 


-2.48 


6.5E-6(0.3) 


40-700 


1.6(0.2) 


3(0.33) 


5,14,23 


991216 


1.02 


317.3(63.4) 


-1.234 


-2.18 


1.9E-4(0.2) 


20-2000 


1.2(0.4) 


4.7(2.3) 


6,13,24,25 


011211 


2.14 


59.2(7.6) 


-0.84 


-2.3 


5.0E-6(0.5) 


30-400 


1.50(0.02) 


3(0.33) 


7,14,26 


020124 


3.2 


110.0(22.0) 


-1 


-2.3 


6.8E-6(0.7) 


30-400 


3.0(0.4) 


3(0.33) 


8,15,27 


020405 


0.69 


192.5(53.8) 





-1.87 


7.4E-5(0.7) 


15-2000 


1.67(0.52) 


3(0.33) 


9,16,28 


020813 


1.25 


211.0(42.0) 


-1.05 


-2.3 


1.0E-4(0.1) 


30-400 


0.43(0.06) 


3(0.33) 


10,17,29 


030328 


1.52 


109.9(21.8) 


-1 


-2.3 


2.6E-5(0.2) 


30-400 


0.8(0.1) 


3(0.33) 


11,17,27 


030329 


0.1685 


67.6(2.6) 


-1.26 


-2.28 


1.1E-4(0.1) 


30-400 


0.5(0.1) 


1(0.11) 


12,18,30 



NOTE. — (a) The spectral parameters fitted by the Band function; (b) the fluence and error observed in the corresponding energy range; (c) the observed break time 
and error of the afterglow light curve; (d) the medium density and error from afterglow fittings, if no available the value of n taken to be 3 ± 0.33 cm ~ 3 ; (e) references 
in order for redshift, spectral data, fj, and n. 

References.— (1) Djorgovski et al. 1998a; (2) Djorgovski et al. 1998b; (3) Hjorfh 1999; (4) Amati et al. 2000; (5) Galama et al. 1999; (6) Vreeswijk et al. 
1999a; (7) Andersen et al. 2000; (8) Hjorfh et al. 2003; (9) Masetti et al. 2002; (10) Price et al. 2002; (11) Rol et al. 2003; (12) Greiner et al. 2003; (13) Jimenez et 
al. 2001; (14) Amati et al. 2002; (15) Barraud et al. 2003; (16) Price et al. 2003; (17) Atteia 2003; (18) Vanderspek et al. 2004; (19) Djorgovski et al. 2001; (20) Frail 
et al. 2003; (21) Kulkarni et al. 1999; (22) Masetti et al. 2000; (23) Bjornsson et al. 2001; (24) Halpern et al. 2000; (25) Panaitescu & Kumar 2002; (26) Takobsson 
et al. 2003; (27) Ghirlanda et al. 2004; (28) Price et al. 2003; (29) Barfh et al. 2003; (30) Berger et al. 2004 
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Fig. 1.— 



Beaming-corrected 7-ray energy versus local-observer peak energy for the GRB sample listed in Table 1. The line is the best fit. 
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FIG. 2. — Hubble diagrams for the GRB sample {filled symbols, for C= 1.12) and the binned SN la data (open symbols) from Riess et al. (2004). The line 
corresponds to a flat cosmology with Qm = 0.27 and h = 0.71. 
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FIG. 3 . — Contours of likelihood in the Qm-^A plane. The cross indicates the best fit, and the dashed line is for a flat universe. 
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FIG. 4. — Contours of likelihood from the GRB sample (dashed lines) in the Qm-w plane. The solid contours consider a prior of S1m = 0.27 ± 0.04. The cross 
indicates the best fit. 



